Reversal-field memory in magnetic hysteresis 
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We report results demonstrating a singularity in the hysteresis of magnetic materials, the reversal- 
field memory effect. This effect creates a nonanalyticity in the magnetization curves at a particular 
point related to the history of the sample. The microscopic origin of the effect is associated with a 
local spin-reversal symmetry of the underlying Hamiltonian. We show that the presence or absence 
of reversal-field memory distinguishes two widely studied models of spin glasses (random magnets). 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

Hysteresis is a widely observed phenomenon. 1 While 
many basic features are well understood 2 *^ a number of 
effects remain unresolved, such as the noise spectrum, ■ 
frequency dependence^ and exchange bias£& in the hys- 
teresis of magnetic recording media. 

Here we report a novel memory effect in magnetic 
systems that emerges when the magnetic field is first 
decreased from its saturation value and then increased 
again from some reversal field Hr. We find that the sys- 
tem exhibits a singularity at the negative of the reversal 
field, —Hr, in the form of a sharp kink in the magnetiza- 
tion. Microscopically a finite number of "symmetric clus- 
ters" are required for this effect to take place. In these 
clusters the central spins flip after all spins on the cluster 
boundary have flipped. In addition, the central spins ex- 
perience zero effective local field so they are symmetric 
with respect to the change of direction of the external 
field. We have observed this "reversal-field memory" ef- 
fect numerically in spin glasses, as well as experimentally 
in magnetic thin films commonly used in recording me- 
dia. 

To better illustrate the reversal-field memory effect 
we use the recently introduced first order reversal curve 
(FORC) diagram approach 9 which captures the distribu- 
tion of characteristic properties of hysteretic systems in 
great detail. 

II. MODEL AND ALGORITHM 

We study the Edwards- Anderson (EA) Ising spin-glass 
HamiltonianiS 

h y^-l.s.s, //»• (i) 

(i,i) 

Here Si = ±1 are Ising spins on a square lattice of size 
N = L x L in two dimensions with periodic boundary con- 
ditions. The exchange couplings Jij arc random nearest- 
neighbor interactions chosen according to a Gaussian dis- 
tribution with zero mean and standard deviation unity, 
and H is the external magnetic field. The complex energy 
landscape of the EA spin glass will give rise to hysteretic 
behavior. 



We simulate the zero temperature dynamics of the EA 
model by changing the external field H in small steps, 
first downward from positive saturation, and then up- 
ward from a reversal field Hr. After each field step, the 
effective local field hi of each spin Si is calculated: 

h i = Y f J ij S j -H. (2) 
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A spin is unstable if hiSi < 0. We then flip a randomly 
chosen unstable spin and update the local fields at neigh- 
boring sites and repeat this procedure until all spins are 
stable. 

III. RESULTS 

Figure ^ shows a typical reversal curve. The area 
around —Hr is enlarged in the inset and shows a sharp 
"kink" . The observation of any such sharp feature in a 
disordered system, especially of finite size and after dis- 
order averaging, is quite remarkable. 

The size of the change in slope can be characterized by 
measuring the slope of the magnetization curves to the 
left and right of —Hr, and comparing the difference with 
the average (see Fig. 12}. The slope changes abruptly by 
as much as 30 % as we pass through H — —Hr. For 
this model, the range of reversal-field values for which 
the kink is present is roughly — 4.5 < Hr < —1.5. 

To develop an initial understanding of the microscopic 
origin of reversal-field memory, we first describe the con- 
ditions under which the effect occurs within the Preisach 
model^i In that approach, a magnetic system is de- 
scribed as a collection of independent two-state (±1) 
switching units ( "hysterons" ) . Unlike Ising spins, which 
always align with their local field, the hysteron's state 
changes from —I to +1 at a field Hf, + H c different from 
the field H\, — H c required to switch the hysteron from +1 
to —1. Here Hi, is the bias field and H c is the coercivity 
of the hysteron. Different systems are then distinguished 
by their distribution p(Hb,H c ) of hysterons of a given 
bias and coercivity, the so-called "Preisach function." 

For "symmetric hysterons," i.e., Hj, = 0, starting from 
a fully 'up' polarized state and decreasing the field to a 
negative Hr switches down all hysterons with H c < \Hr\. 
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FIG. 1: Reversal-field memory effect: starting from positive 
saturation, the field is decreased to a negative reversal field 
Hr, then increased again. The inset shows the kink at —Hr. 
Data are from a two-dimensional EA Ising spin glass with 
100 2 spins. The dotted line represents the major hysteresis 
loop, the dashed lines mark ±Hr, respectively. In all figures 
the error bars are smaller than the symbols and the data are 
averaged over 1000 disorder realizations. In addition, H and 
M are in units of the interaction bond standard deviation, 
which in this work is set to unity. 




FIG. 2: The difference (triangles) and average (circles) of the 
left and right derivatives on either side of —Hr for the EA 
Ising spin glass. 



Reversing the direction of the sweep and increasing the 
field from Hr to — Hr switches back every switched hys- 
teron. Thus at H = —Hr saturation is reached, creating 
a kink in the magnetization. Symmetric hysterons there- 
fore give rise to reversal-field memory, but only if their 
number is macroscopic. In this case p(Hb, H c ) must have 
a Dirac delta singularity at Hf, — and H c = \Hr\. 
If the kink is observable in a range of Hr values, as in 
Fig. then the singularities of the Preisach function form 
a ridge along the Hi, = axis for that range of H c values. 
The spins in the EA model are not independent. For 
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FIG. 3: Overlap function q of the spins flipping at Hr and 
at H > Hr, for Hr — —2.28. The insets show data for 
Hr — —0.40 and Hr — —5.60, respectively. The arrow in the 
inset marks H — 0.40. 



them to behave as symmetric hysterons, they must pos- 
sess a local spin-reversal symmetry. By local spin- 
reversal symmetry we mean that the local field hi, felt 
by Si [Eq. ©], is perfectly reversed if the external field 
H is reversed and all spins coupled to Si are reversed as 
well. In the EA Hamiltonian a cluster has this symmetry 
when all neighbors of Si align with the local field before 
Si does, both for decreasing and increasing fields, i.e., Si 
flips down only after all of its neighbors are already neg- 
ative, and during the reverse sweep Si switches up only 
after all the neighbors. 

The importance of the symmetric clusters can be ap- 
preciated by noticing that every spin of the EA model 
seems to have local spin-reversal symmetry. However, 
note that the spin configurations in general depend on 
the history of the sample in the glassy phase. Therefore, 
at —Hr the neighbors of most spins do not necessarily 
point in a direction opposite of their direction at Hr, 
thus most EA spins do not belong to symmetric clus- 
ters. Hence the model Hamiltonian possessing a local 
spin-reversal symmetry is a necessary but not sufficient 
condition of having symmetric clusters. In order to test 
this conclusion we studied the random field Ising model 
(RFIM)iSiiii In this model the couplings Jy are chosen 
to be uniform, and the disorder is introduced through 
random local fields hi . Direct inspection reveals that the 
RFIM does not possess a local spin-reversal symmetry, 
nor did we find a magnetization kink in our simulations. 

To further illustrate the reversal-field memory effect, 
we define an overlap function q between the spins which 
flip at Hr and the spins which flip at H > Hr: 

= \ Y\ Si ( HR + *) - Si ( H *y x 

i 

[Si(H + 6) - Si(H)) . (3) 
Here 8 is the field step. In Fig. [3] we show the overlap 
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FIG. 4: FORC diagram of the EA Ising spin glass. Note the 
ridge along the H c axis. 
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FIG. 5: Experimental FORC diagram of a Kodak sample. 
Note the similarity to the FORC diagram of the EA Ising 
spin glass shown in Fig. 0] 

versus H > Hr for Hr — —2.28. Clearly the spins which 
flip at Hr are highly correlated with the spins which flip 
at —Hr. This in turn shows that we have a macroscopic 
number of symmetric clusters. The insets show a much 
weaker effect at Hr = —0.40 and Hr = —5.60, values 
outside the ridge of Fig. [5] 

The reversal-field memory can be characterized in 
greater detail by adapting a new tool developed for ana- 
lyzing experimental data for hysteretic systems^ In this 
method, a family of FORCs is generated by varying Hr 
in small steps. Let us denote by M(H,Hr) the mag- 
netization as the function of the applied and reversal 
fields. Computing the mixed second derivative with re- 
spect to H and the reversal field Hr yields the "FORC 
distribution" p(H c ,Hf,), where H c = (H — Hr)/2 repre- 
sents the local coercivity and H^ = (H + Hr)/2 the bias. 
FORC distributions are more general than Preisach dis- 
tributions since they are model independent and allow 
for asymmetries and negative regions. 



Figure 0] shows the FORC diagram of the EA model. 
The ridge along the H c axis corresponds to the peak of 
Fig. [21 and is therefore a striking representation of the 
reversal-field memory effect. In contrast, the RFIM ex- 
hibits a vertical structure, even though both models have 
quite similar major hysteresis loops. This sensitivity of 
the FORC method makes it uniquely suited for analyzing 
details of hysteretic systems. 

To demonstrate the existence of the reversal-field mem- 
ory in magnetic media, in Fig. [S] we show an experimen- 
tally obtained FORC diagram of well-dispersed nonin- 
teracting single-domain magnetic Co-7-Fe203 particles. 
Typically FORC diagrams of various systems exhibit ex- 
tensive differences (see for example Ref. IT2T) . There is 
a striking similarity between the experimentally and nu- 
merically determined FORC diagrams: both exhibit a 
narrow ridge along the H c axis, the ridge "melting" with 
decreasing H c and narrowing with increasing H c . Sim- 
ulations on more realistic systems with dipolar interac- 
tions, such as the ones found in the experimental sample, 
also show evidence of a ridge and therefore reversal-field 
memory. 

IV. CONCLUSIONS 

In conclusion, we report a reversal-field memory ef- 
fect in disordered magnetic systems that manifests it- 
self as a sharp kink in first order reversal curves, and 
as a sharp ridge on the zero bias axis of FORC dia- 
grams. The effect has been observed numerically in the 
Edwards-Anderson Ising spin glass and experimentally 
in well-dispersed recording media samples. The origin of 
this effect is tied to a local spin-reversal symmetry, and 
has its microscopic origin in the development of a macro- 
scopic density of symmetric clusters. Our studies also 
establish that the FORC method is a powerful diagnos- 
tic tool for characterizing magnetic materials. 
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